The Structure of the Big Bang from 
Higher-Dimensional Embeddings 
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We give relations for the embedding of spatially-flat Friedmann-Robertson- Walker cosmological 
models of Einstein's theory in flat manifolds of the type used in Kaluza-Klein theory. We present 
embedding diagrams that depict different 4D universes as hypersurfaces in a higher dimensional flat 
manifold. The morphology of the hypersurfaces is found to depend on the equation of state of the 
matter. The hypersurfaces possess a line-like curvature singularity infinitesimally close to the t — + 
3-surface, where t is the time expired since the big bang. The family of timelike comoving geodesies 
on any given hypersurface is found to have a caustic on the singular line, which we conclude is the 
5D position of the point-like big bang. 
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I. INTRODUCTION 



In four-dimensional general relativity as formulated by 
Einstein, the singularity theorems tell us that the big 
bang was a unique birth event for a universe whose mat- 
ter is of the conventional sort and obeys the energy condi- 
tions H], D . The early stages of the standard Friedmann- 
Robertson- Walker (FRW) models may be modified by in- 
flation, basically because then there is a period when the 
energy conditions are violated || Q . The associated pro- 
cess of mass generation can either be described at the par- 
ticle level by quantum field theory or at the macroscopic 
level by classical theory []|, ||, [?], [| ^, [rJ . The early stages 
of standard cosmology may also be modified by the ef- 
fects of quantum gravity and quantum tunneling |l2| . 
In this regard, the appropriate formalism in 4D for early 
times can be investigated by considering the problem in 
5D p3[ . Also, the status of standard cosmological models 
for late times can be clarified by considering how they are 
embedded in 5D |l4j]. There has recently been renewed 
interest in 5D manifolds of the type used originally in 
Kaluza-Klein theory 



15 . This is partly because of 



the resurrection of Campbell's theorem, which ensures a 
perfect local embedding of a 4D Einstein space in a 5D 
Kaluza-Klein space |r| [l7| and partly because of 
the implications of higher-dimensional approaches such 
as superstrings and supergravity fl^, p0, Ell |2^] as well 
as string and membrane theory [E3[ E4, Recently, 
several authors have revisited the problem of embedding 
4D FRW metrics in 5D Minkowski space with an eye to- 
wards determining the geometrical structure of the big 
bang (2(| For the standard spatially-flat (k = 0) 

4D FRW models the appropriate class of 5D models was 
given in a notable paper by Ponce de Leon [E8j. His 
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models are separable in time t, space a, and an extra co- 
ordinate I, and reduce to the FRW ones on I = constant 
hypersurfaces (henceforth denoted by Ylg). Interesting 
work has been done on those and related models, both 
from the astrophysical side ||^, ||(| [H| and the mathe- 
matical side ]3^, [33], |34|. One of the Ponce de Leon 5D 
solutions when unembedded is the 4D Einstein-de Sit- 
ter solution, which is in tolerable agreement with obser- 
vations; but surprisingly this and other Ponce de Leon 
solutions have been shown by fast algebraic computer 
package |35| to be flat in 5D. This may appear unusual, 
but recent work has made it clear that a flat and ap- 
parently em pty 5D manifold may contain a curved, 4D 
manifold [hi], [3^, |33j. ^g|, [37|. However, while much 
work has been done on the 5D Ponce de Leon solutions, 
no systematic account has been given of how they embed 
the 4D big bang. 

The present paper addresses this issue. The plan of 
the paper is as follows: Section II introduces the Ponce 
de Leon class of one-parameter 5D cosmological solutions 
to the vacuum Einstein field equations. hypersurfaces 
are shown to be geometrically identical to spatially flat 
(k = 0) FRW cosmological models. The type of matter 
in these complementary 4D models is determined by the 
value of the parameter a, and includes both inflation- 
ary and non-inflationary cases. The 5D manifold is then 
shown to be identical to Minkowski space M5 via a coor- 
dinate transformation. This coordinate transformation 
allows us to plot our universe as a hypersurface in a 3D 
Cartesian space once angular variables are s uppre ssed, 
which is the subject of section III. Subsection III A deals 



with analytically determining where comoving observers 
on Y.£ emanate from in M5 as t — > + , which corre- 
sponds to a singularity in the 4-geometry. Subsection 
III B demonstrates that the local geometry of £f around 



regular points is elliptical for non-inflationary cases and 
hyperbolic for inflationary ones. Subsection III C estab- 
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lishes some results concerning the global properties of the 
foliation of M5, the main conclusion being that the 
cosmological coordinates used in the Ponce de Leon met- 
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ricdo not cover all of 5D Minkowski space. Subsection 
[IIP presents computer plots of T, e hypersurfaces for a 
wide varie ty of cases. These plots are discussed in subsec- 



tion [HE . Also in this subsection, a detailed analysis of 
the singular points on is presented. We find that there 
is a line-like curvature singularity on E^, but that the big 
bang corresponds to only one point on that line. Section 
IV discusses special cases of the Ponce de Leon class not 
mentioned in other sections and the corresponding 4D 
universes. Section [v| contains our concluding remarks. 
Conventions. Throughout the paper, the 5D metric 

signature is taken to be (H ), while the choice 

of 4D metric signature is (H ). The spacetime co- 
ordinates are labeled a: = t, x 1 = (r,8,<fi). The extra 
coordinate is x A — £. The range of tensor indices is as 
follows: A,B,... = 0-4; a,/3,... = 0-3; and a,b = 0, 1. 

indicates the 5D covariant derivative. We denote 
the metric of fiat 3-space by da 2 — dr 2 + r 2 dil 2 with 



dfl 2 = d6 2 + sin 2 
where c = G = 1. 



Finally, we use geometric units 



II. PONCE DE LEON COSMOLOGIES 

The standard class of Ponce de Leon cosmological solu- 
tions is flat in ordinary 3D space, curved in 4D spacetime 
and flat in 5D g g |o|, ^ §2[ || |§. It has a 5D 
line element 

dS 2 = I 2 dt 2 - * 2 /^2/(i-°0 do 2 - a 2 (l - a)- 2 t 2 dl 2 . 

(1) 

This metric is an exact solution of the 5D field equations 
in apparent vacuum, which in terms of the Ricci tensor 
are Rab = 0. In what follows, we will examine the E^ 
hypersurfaces of (|l|) in detail. The induced metric h a p 
on these hypersurfaces is given by 

ds 2 = h aP dx a dx S) = l 2 dt 2 - t 2 ' a i 2 l^- a Ua 2 , (2) 

which we can, of course, identify with the line element of 
a spatially-flat (k = 0) FRW cosmology with scale factor 



fl(r) = *(2«-l)/a(l-a) T l/a 



(3) 



where r = it corresponds to the FRW clock time. That 
is, the geometry of each of the E^ 4-surfaces is identical 
to geometry of the 4D spacetime models commonly used 
to describe our universe. 

We can use the form of the induced metric (Q) to cal- 
culate the 4D Einstein tensor G a p on each of the E^ 
surfaces. Alternatively, we can use standard techniques 
to break the 5D field equations Rab — into one wave 
equation, four conservation equations and ten equations 
for the components of G a p j|; [l6|. The latter are 
given in general by 



= — ^ r TT^T 
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where (744 = — $ 2 , and a star (*) denotes d/dt. We now 
ask the question: what kind of 4D matter would give 
rise to such an Einstein tensor? To answer, we must 
associate G a p with a stress-energy tensor T a p via the 
Einstein equations G a /3 — &irT a p. For the Ponce de Leon 
solutions, this "induced" stress-energy tensor is of the 
perfect-fluid type with the density and pressure of the 
"induced" matter given by: 



8irp 



87175 



2a -3 



(5) 



It should be stressed that this density and pressure refers 
to a 4D matter distribution that gives rise to a curved 4D 
manifold, via the Einstein equation, geometrically iden- 
tical to the hypersurfaces. We have not inserted any 
matter into the 5D manifold by hand, as is commonly 
done in the brane- world/thin-shell formalism where the 
5D stress-energy is a discontinuous, distribution-valued 
tensor field. It is for this reason that we call the matter 
described by T a p induced; rather than being inserted into 
the Einstein equations as an external source for the gravi- 
tational field, the matter distribution is fixed by the mor- 
phology of the surfaces and the higher-dimensional 
geometry. 

The equation of state of the induced matter is 



p= (2a/3-l)p, 



(6) 



which is determined by the dimensionless parameter a in 
the metric (0). For a = 3/2, 87179 = 4/3r 2 and p = 0, 
which describes the late (dust) universe. For a = 2, 
87173 = 3/4r 2 and 8717? = l/4r 2 so p — p/3, which de- 
scribes the early (radiation) universe. That these quan- 
tities are the same as they are in the 4D FRW models 
should not surprise us, as this could be inferred from 
Campbell's theorem [l4|, [r|. However, the Ponce de 
Leon class of solutions also embeds some more exotic 
cosmological scenarios. To see this, let us form the grav- 
itational density of the induced matter, 



p g = p + ip 



3(a - 1) 



(7) 



The strong energy condition demands that p g > 0, which 
is only satisfied for a > 1. For cases which violate this 
condition a € (0, 1), the Raychaudhuri equation says 
that geodesies paths on accelerate away from one an- 
other, as is expected in inflationary scenarios. This is 
easily verified by examining the congruence of comoving 
da = paths in the 4-geometry described by (|J). Just 
like in conventional FRW cosmology, the proper distance 
between adjacent paths will increase at an accelerating 
pace if the deceleration parameter, 



. . d 2 a ( da 

q{T) -- a d^{d^ 



1. 



(8) 



is negative. We therefore conclude that models with 
a S (0, 1) describe inflationary situations where the cos- 
mological fluid has repulsive properties, while models 
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with a > 1 have ordinary gravitating matter. Note that 
we exclude models with a < because they imply a con- 
tracting universe. 

As a quick aside, we note that the Ponce de Leon met- 
rics (^) have a constrained equation of state (i.e., the type 
of matter corresponding to the curvature of the E^ hyper- 
surfaces does not change in time) . If one wanted to study 
a more realistic model of our universe, the radiation- 
matter transition could be generated by joining metrics 
with a — 2 and a = 3/2 across a t = constant hypersur- 
face that represents the surface of last scattering. Such a 
calculation is, however, beyond the scope of this paper. 

It was discovered by computer that (|l|) is actually flat 
in 5D {H |0|, II || H, || |H. To prove this alge- 
braically is non-trivial, but we can change coordinates 
from those in (Fn) to 



T(t,r,£) = 



" 'i + r_wv^v(i-») 



2 \ a 
2 



1 - 2a 
R(t,r,i) = rt x l a i x l^- a \ 
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Then Oil) becomes 



dS 2 = dT 2 



(dR 2 + R 2 dn 2 



dL 2 



(9a) 
(9b) 

(9c) 
(10) 



which is 5D Minkowski M5. In order to contrast with 
the standard coordinates x A — (i, r, 8, cj), £), we will label 
the Minkowski coordinates as y A = (T,R,8,(j),L). The 
importance of ([)]) is that it allows us to visualize the 
geometric structure of our universe when viewed from 
5D flat space by plotting the surfaces defined by t = 
constant. As seen in equation (^|) above, these hyper- 
surfaces share the same intrinsic geometry as standard 
4D FRW k = cosmologies, which means that they 
are equivalent to such models in a mathematical sense 
(this is the paradigm commonly adopted when one tries 
to embed known solutions of the Einstein equation in 
higher-dimensional manifolds 26, 27|). We can therefore 
learn a great deal about the local and global topological 
properties of 4D FRW models, as well as the geometric 
structure of the "physical" big bang, from studying the 
E^ hypersurfaces. These and other issues are considered 
in the following sections. 



III. THE STRUCTURE OF THE EMBEDDED 
UNIVERSES 

In order to analyze the geometrical properties of the 
Ef hypersurfaces, we note that the T, R, L coordinates 
of ( |l0| ) are orthogonal and so can be used as Cartesian 



axes. Then, the universe will have a shape defined by 
T = T(t,r,£ ), R = R{t 7 r,£ ), L = L(t,r,i ) on some 
hype rs urface I = £q. The following three subsections 
ilLA| - |IIIC| explore various characteristics of the E^ hy- 
persurfaces analytically. The goal is to get a feel for the 
geometry of the E^ hypersurfaces before computer plots 
are presented in subsection III D , In that subsection, we 



plot hypersurfaces for several different cases, suppressing 
the and (j) coordinates so that our universe appears as a 
2-surface embedded in Euclidean 3-space E3. Subsection 
HIE deals with the location and five-dimensional shape 



of the big-bang singularity, and is in some sense a con- 
tinuation of the discussion of singular points started in 
subsection III A. 



A. Properties of Singular Points on Ef 

In this subsection, we will begin to investigate singular- 
ities on the Ef hypersurfaces, and introduce the problem 
of locating where these singularities occur in M5. Our 
goal is to get a feel for the shape of the E^ surfaces us- 
ing analytical method s befo re resorting to the computer 
plotting of subsection HID . 

A cursory inspection of equation (||) reveals what ap- 
pears to be a singularity in the 4-geometry at t = 0, 
which is commonly associated with the big bang. We 
can calculate the 4D Kretschmann scalar (or curvature 
invariant) of the £ = £q hypersurfaces: 



_ pQ/375 



K = R 



12 (a 2 -2a + 2) 



(11) 



We note that K diverges for all a at t — 0. This implies 
that the anomaly at t — is a genuine curvature singu- 
larity in the 4D manifold represented by However, 
we know that the 5D manifold ([!]) has Rabc d — and is 
therefore devoid of any singularities. Therefore, the 4D 
big bang results from the nature of the E^ hypersurfaces, 
or equivalently the choice of 5D coordinates. 

The question naturally arises of where the big bang at 
t = appears in M5 . One way to approach the problem is 
to take the t — > + limit of the coordinate transformation 
(^) while holding r constant. It transpires that there are 
two different physically interesting cases: 



lim T(t,r,£) 



lim R(t,r,£) 
t->o+ 



lim L(t,r,£) 



£a/(l-a) 

2(1 -2a) [ t 

Jo, a <E (5,00) 
1 -00, a S (0, |) 

0, a > 

£a/(l-a) 



(l-2a)/a 



(12a) 
(12b) 



2(1 - 2a) 



t 



(l-2a)/c 



0, a € (i, 00) 
+00, ae(0,i) 



(12c) 
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In other words, the congruence of (r, 9, <fi) — constant 
curves on E^ (henceforth denoted by j r ) converges to a 
point as t — > + for all a > 0. For a 6 (|, oo) the caustic 
is at T = i? = L = 0; for a € (0, ^), the caustic is at null 
infinity as approached by following the ray T+L = R = 
into the past. Naively, this would seem to suggest that 
the the big bang is a point-like event in 5D. However, 
caution is warranted because the limiting procedure de- 
scribed above is not unique. In particular, we need not 
approach t — in a manner that leaves r = constant. For 
example, suppose we follow the path r(t) oc t~ x l 2a into 
the past. If a £ (^, oo), then we will end up at the 5D 
point T = —L = constant ^ 0, R = 0, which is different 
from the caustic in the "f r congruence. So, it is not cor- 
rect to uniquely associate t = with the position of the 
7 r caustic. In general, there are many points on E^ eligi- 
ble for that distinction. On the other hand, it should be 
stressed that if we approach the initial singularity along 
any path where r(t) tends to a finite constant as t — » + , 
then we will end up at one of the 7 r caustics (which one, 
depends on the value of a). The issue of the precise lo- 
cation of the big bang in 5D is compl icated , so we will 



where ^gf b is the Euclidean 2-metric on S and ( 2 >Rf, 



1212 



defer a full discussion until subsection [IIIE 

Regardless of where the t — surface lies in M5, 
the calculations presented in this section demonstrate 
that the point T=i? = I = 0isa special point 
for a £ (i,oo), and that a point at null infinity T — 
—L — ► —00, R = is a special point for a £ (0, ^). The 
computer-generated figures presented in subsection [II D 
will bear out this analytical conclusion. 



B. Properties of Regular Points on 

Having already given a partial discussion of the singu- 
lar points on the E^ hypersurfaces, we would now like to 
concentrate on analytically determining what E^ looks 
like in the neighborhood of regular points. As mentioned 
above, our plots of E^ will look like 2-surfaces embedded 
in E3. It is therefore useful to recall the notion of the 
Gaussian curvature G E of 2-surfaces S £ E3 from stan- 
dard differential geometry J38|. Consider a point P on 
a surface S embedded in E3. In a small neighborhood 
around P, the surface 5* can be modeled as a paraboloid 
centered at P (the so-called "osculating paraboloid"). 
The shape of this paraboloid is given by the sign of the 
Gaussian curvature, which in turn is given by the prod- 
uct of the principle curvatures of S at P. If G E > 0, 
the paraboloid is elliptical and the surface has a convex 
or concave shape about P. If G E < 0, the paraboloid is 
hyperbolic and the local shape of S is that of the familiar 
"saddle-surface". If G E = 0, one or both of the princi- 
ple curvatures are zero implying a cylindrical or planar 
shape. In terms of intrinsic geometrical quantities, the 
Gaussian curvature is given by 



G = 



(2) <oi 
detp) ff E b ] 



(13) 



is the single independent component of the associated 2D 
Riemann-Christoffel tensor. For the Euclidean version of 
the Ef hypersurfaces, we can find ^gf b from the positive- 
definite line element 



(2) 



gl b dx a dx b = (dT 2 + dR 2 + dL 2 



d£=Q 



(14) 



Using this, we can calculate the Gaussian curvature ( fl3| ) 
of the E^ surfaces when embedded in Euclidean space: 

jj _ Aa 6 (a - l)i2(a-2)/a £ -2(l+a)/(l- Q ) 
G — : — : : ; — : , (15a) 



where 



[x 2 — 2r 2 x + {r 2 + a 2 ) 2 



X = a 2 i 2(Q_1)/Q £ 2Q/(Q_1) . 



(15b) 



We note that the denominator of G E is positive definite 
for all values of t, r, £ and a. Therefore, G E > for 
a £ (l,oo) and G E < for a £ (0,1) for all (t,£) > 0. 
In other words, when viewed in E3, all the non-singular 
points on E^ are elliptical for non-inflationary models and 
hyperbolic for inflationary ones. 

We mention in passing that it is also possible to calcu- 
late the Gaussian curvature G M of the E^ surfaces when 
they are embedded in M3. The 2-metric ^g^ b is obtained 
by setting d6 = d(f> = dt = in (Eh. This yields 



WR 
1 i 



0101 



1 



det[( 2 ) ff M] a 2 ^ 2 ' 



(16) 



which, like G E , is positive for a £ (1, 00) and negative 
for a £ (0, 1). Also, we note that G M diverges as t — > + 
for all a, which is consistent with the behaviour of the 
4D Kretschmann scalar (O ) calculated above. The Ricci 
tensor in 2D is given by^i?jf b = (2) gf fc G M , so that the 
strong energy condition reads ^ R^f b u a u b = G M > 
where u a is some timelike 2-velocity. Hence, if G M < 
then the strong energy condition is violated in 2D. This 
implies that the 2D gravitational density of the matter is 
less than zero, which in turn means that the relative ve- 
locity between neighboring timelike geodesies is increas- 
ing in magnitude by Raychaudhuri's equation. This is 
consistent with the finding that q(t) < and p + 3p < 
for a £ (0, 1) discussed in subsection II. 

In conclusion, we have found out that the local geom- 
etry around regular points on E^ - for a given value of a 
- is determined by whether or not the model represents 
a universe with ordinary or inflationary matter. In the 
former case the regular points are elliptical, in the latter 
case they are hyperbolic. 



C. Global Properties of 

In the previous two subsections, we have discussed the 
properties of singular and regular points on E^. We now 
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turn our attention to the global properties of the hyper- 
surfaces, with an eye to determining which portion of M 5 
they inhabit. In other words, we want to know what part 
of the M5 manifold is covered by the x A = (t,r,8,cj),£) 
coordinates. 

The coordinate transformation (|^) implies 



(T + L)/R= a/r. 



(17) 



Since R/r is positive for (t,£) > and we assume that 
a > 0, this implies that the Tig hypersurfaces are re- 
stricted to the half of the manifold defined by (T+L) > 0. 
Next, the transformation also implies that 



T - R 2 — L = 



2a - 1' 



(18) 



Hence, if a G (i, 00) then the Tg surfaces must lie in the 
region satisfying T 2 — R 2 — L 2 > 0, which is the volume 
inside the light cone originating from T = R = L = 0. If 
a G (0, 5), the surfaces must lie in the region outside the 
light cone, defined by T 2 - R 2 - L 2 < 0. Also, the Tg 
surfaces approach T 2 — R 2 — L 2 = as I — > 0, implying 
that the direction of increasing £ is always away from 
the light cone. Putting these facts together, we see that 
if a G (|,oo), then the coordinates utilized in (Q) only 
cover the portion of M5 inside the light cone centered at 
the origin with T > 0. If a G (0, h), then the coordinates 
in |l]) cover the region of M5 satisfying (T + L) > and 
exterior to the light cone centered at the origin. Finally, 
equations ( pb| ) , ( p"7j ) and ( jig ) can be rearranged to yield 



= UV 



a 2(l-Q)^2(l-2Q)/(l-Q) [/ 2a 



(2a -1) 



R\ 



(19) 



where we have introduced the advanced and retarded co- 
ordinates 



U = T + L, V = T-L, 



(20) 



This defines the Tg hypersurfaces entirely in terms of the 
y A coordinates and I. This is a useful relationship which 
we will need below. 

To conclude, we have discovered that the Tg hypersur- 
faces are constrained to lie in one half of M 5 for all a > 0. 
Furthermore, they must lie within the light cone centered 
at the origin for a G (|, 00) and outside it for a G (0, -|). 
We have summarized all that we have learned about the 
effects of a on the Tg hypersurfaces in Table |. 



D. Visualization of the hypersurfaces 

We are now in a position to discuss Figures |l| - [| in 
detail. The plots show the constant t and r coordinate 
lines associated with a given Tg hypersurface (denoted 
by 7( and j r respectively) as viewed from the (T, R, L) 
space. As an aid to visualization, we choose to let r 
and R range over both positive and negative numbers, 




FIG. 1: £f hypersurface with a = 3/2 and Iq = 1. We 
take t G [10" 400 ,2] and r G [-10,10], with At « 0.01 and 
Ar w 0.003. 



so the hypersurfaces are symmetric about the R = 
plane (if it is desired to have r > or R > one of the 
symmetric halves can be deleted). As the models evolve 
in i-time they will generally grow in the i?-direction, so 
that the hypersurfaces appear to "open-up" in the 
direction of increasing T. The 74 isochrones are seen to 
wrap around the hypersurfaces such that they cross the 
R = plane perpendicularly. The j r lines are orthogonal 
(in a Euclidean sense) to 74 lines at R — and run along 
the length of the surface. 

When using a computer to plot the surfaces, it is nec- 
essary to choose finite ranges of t and r. The ranges used 
in each of the Figures are given in the captions, along 
with the increments in the time or radial coordinates be- 
tween adjacent coordinate lines (denoted by At and Ar 
respectively). For models with a G (|,oo), the t = 
"line" is really a point for finite values of r. Therefore, 
the first 7t line in such plots has been chosen to be t = S, 
where S is the smallest positive number allowed by ma- 
chine precision. For such models, the lower edge of the 
Tig hypersurfaces corresponds to the t = 8 line. For plots 
with a G (0, 5) we have chosen to restrict t and r to a 
fairly narrow ranges in order to facilitate visualization. 
As a result of this, the lower edges of Sf in such plots 
are the 7 r curves with the largest value of |r|. How- 
ever, we have confirmed that if the range of r is larger 
and min(t) = 5, then the comoving trajectories with the 
largest values of |r| tend to bunch up along the t = S 
line. That is, if the range of r is unrestricted then the 
lower edge of the T,g hypersurfaces will correspond to the 
7t line closest to the big bang, irrespective of the value 
of a. 

We have carried out an extensive analysis of the mor- 
phology of the Tig hypersurfaces in the (a,£o) parameter 
space, and present five informative cases which are illus- 
trated in the Figures: 
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a interval 


sign of q(r) 


osculating paraboloid 


caustic location 


T 2 -R 2 - L 2 


(o,|) 




hyperbolic 


null infinity 




(1,1) 




hyperbolic 


origin 


+ 


(1,00) 


+ 


elliptical 


origin 


+ 



TABLE I: Dynamical and topological properties of cosmologies embedded on hyp ersurf aces. The "caustic location" refers 
to the point that the ~/ r lines approach as t — ► + . 




Model 1: a = 3/2, £q — 1. This is the standard matter- 
dominated model for the late universe (see above). 
The equation of state is that of dust and the scale 
factor evolves as i 2 / 3 . As expected, the osculating 
paraboloid has a elliptical geometry for every non- 
singular point on Y,g. The congruence of j r curves 
has a caustic at T = R = L = 0. Notice that 
the rate of expansion of j r lines seems to slow as 
time progresses and that the surface lies within the 
light cone centered at the origin. This plot is very 
similar to the embedding diagram shown in Figure 
5 of Lynd cn-Bc ll et al. 26 (see the discussion in 
subsection III E| ) . 



Model 2: a = 2, £q = 1. This is the standard radiation- 
dominated model of the early universe (see above) . 
The scale factor evolves as i 1 / 2 , the caustic in the 
7 r congruence is located at the origin, and Hi is 
qualitatively similar to Model 1. 

Model 3: a = 2/3, Iq = 1. A moderately inflationary 
model with (p+3p) < 0, a scale factor cx i 3 / 2 , and a 
caustic in the j r congruence at the origin. The ge- 
ometry of the osculating paraboloid is hyperbolic, 
7 r lines accelerate away from each other, and the 
entire surface lies within the light cone centered at 
the origin. 



Model 4: a = 1/30, 



1. A very inflationary model 



with (p + 3p) < 0, a scale factor cx t 30 , and the 
"f r caustic at past null infinity. The osculating 
paraboloid is hyperbolic with comoving paths fly- 
ing apart. The surface lies outside the light cone 
centered at the origin. 

Model 5: a = 1/3, i = 20, 40, 60. A set of inflationary 
models with (p + 3p) < 0, a scale factor cx i 3 , and 
the 7 r caustic at past null infinity. As expected, all 
the surfaces lie outside the central light cone. The 
space between hypersurfaces and the null cone 
increases with increasing I. 

We note that the last case is based on classical field 
theory but bears a striking resemblance to computer- 
generated models of stochastic inflation |3^] based on 
quantum field theory. 



E. Comments Regarding the Figures and the 
Geometrical Nature of the Big Bang 

One of the most interesting features of Figures [l] - || 
is that the constant-T cross sections 7t of the hyper- 
surfaces are very nearly closed. We will now demonstrate 
that the finite amount of space between the lower edges 
of the surfaces depicted in the Figures is a result of the 
requirement that t > (5, which is imposed by limitations 



FIG. 4: T,i hypersurface with a = 1/30 and £ = 1. We 
take t G [0.9,1.1] and r G [-0.1,0.1], with At w 0.01 and 
Ar » 0.003. 



FIG. 5: E £ hypersurfaces with a = 1/3 and 4 = 20,40,60. 
We take t G [10" 300 ,3] and r G [-2,2], with At « 0.05 and 
Ar « 0.07. 



in machine arithmetic. If r is eliminated from equations 
(pa[) and (9b), we find that the projection of 7t lines onto 
the T_R-plane is parabolic. For early times, we obtain 



T[2ai 1/c ^ 1/(1 ~ Q) ] = R 2 



1 - 2a 



t < 1. 



(21) 



If we hold T constant and let t — > + then i? — > 0, which 
establishes that the 7t cross sections are asymptotically 
closed. Furthermore, since the lower edges of the surfaces 
correspond to the minimum value of t, the 74 isochrones 
must approach the R = plane as t — > + . By exam- 
ining equations (|J|) and (|l9"|), it is easy to see that E^ 
itself must approach the T + L = plane as i + . 
Therefore, the 7t lines must approach some subset of the 
T + L = R = null line AT as t — > + , which we denote 
by the intersection of A/" and E^: TVnE^. In other words, 
the finite gaps shown in the Figures are an artifact from 
numerical calculations. 

This raises a puzzling and subtle issue. The "ft curves 
approach a line in 5D at t — > + , yet the coordinate trans- 
formation (||) suggests that the t — "line" is really a 
point in 5D for any finite value of r. Is the big bang a 
point-like or line-like event in 5D? Several authors have 
appreciated this problem in the past. Lyndcn-Bcll et al. 
present an embedding diagram of a spatially-flat FRW 
model with p = 0, which is equivalent to our Model 1 
shown in Figure [l] [^6| . Their plot is qualitatively similar 
to ours, and they identify the half of M satisfying T > 
with the 4D big bang. This is contrary to the conclu- 
sion of Rindler, who demonstrates that open FRW mod- 
els may be completely foliated by spacelike 3-surfaces of 
finite volume |27|| . The volume of these hypersurfaces 
tends to zero as they approach the initial singularity, 
which leads to the idea that the big bang is a point-like 
event. 

We will now demonstrate that both of these views are 
in some sense mathematically correct; the location of the 



big bang depends on one's definition of the initial singu- 
larity. One such definition might be the locus of points 
defined by the positions of the fundamental comoving ob- 
servers in the t — > + limit. This is in the spirit of the def- 
inition of a singular spacetime as a manifold containing 
one or more incomplete timelike geodesies; i.e. timelike 
geodesies which are inextendible in at least one direc- 
tion, or have a finite affine length 0). The incomplete 
geodesies on are clearly the comoving trajectories of 
the 7 r congruence, which has a fundamental caustic as 
t — » + . For a G (5,00), it is obvious from Figures [j] 
- U that it is impossible to extend these geodesies past 
the caustic at T = R = L = 0. The 7 r geodesies are 
also incomplete for a S (0, 5), although this is hard to 
see from Figures || and |H| because the caustic is at null 
infinity. However, we note that the proper time interval 
along a comoving path between the initial caustic and 
any point on T,£ is finite, which establishes that the 7 r 
curves are incomplete. In both cases, the comoving paths 
radiate from a point in 5D. By identifying the location 
where the 7 r congruences "begin" with the initial singu- 
larity, we conclude that the big bang is a point-like event 
in higher dimensions. 

However, the definition of the big bang need not be 
tied to the properties of inextendible timelike geodesies. 
We could instead elect to define the location of the big 
bang as the locus of points on where curvature scalars 
diverge. This definition is in the spirit of using quanti- 
ties like the Kretschmann scalar to distinguish between 
coordinate and genuine singularities in the Schwarzschild 
spacetime (among others). For the present case, we note 
that the Copernican principle built into the FRW so- 
lutions demands that curvature scalars on must be 
constant along 74 isochrones. This is evidenced by the 
Kretschmann scalar (|ll]) and the Minkowski Gaussian 
curvature (|l^), both of which demonstrate no depen- 
dence on r. As t — > + , both of these scalars diverge 
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and the 74 isochrones approach TV PI E^. This would seem 
to suggest the line TV" n is the location of a curvature 
scalar singularity on By the definition mentioned at 
the head of this paragraph, this implies the big bang is a 
line-like event in higher dimensions. 

There are two possible objections to this conclusion: 
Firstly, the idea that E^ is singular along NnS; is some- 
what counterintuitive, because the hypersurfaces shown 
in the Figures appear to smoothly approach TV as t — > + . 
Secondly, we observe that the (t,r) coordinates do not 
actually cover TVn E^, which makes the ^-dependent ar- 
guments of the last paragraph suspect. This is easily 
seen by noting that the induced metric (^|) is singular at 
t = and T + I = f = Oon/VnEf from equations 
( pb| ) and (0). The situation is analogous to the failure 
of spherical coordinates to cover the z-axis in E3. For 
these two reasons, it would be a good idea to confirm 
the presence of a curvature singularity along TV n E^ in a 
manner independent of the (t, r) coordinates. 

This can be accomplished using the well-known extrin- 
sic curvature formalism. Restoring the 9 and cj> coordi- 
nates, let us now work in the 5D y A — (U, R, 8, </>, V) 
system with U and V as defined in (|20|). In this coordi- 
nate system, the line TVfl E^ is approached by taking the 
U — » + limit (recall that U > for all of the hypersur- 
faces). We can define E| in this space by 



= $t(U,R, V), 



(22) 



where $f (U, R, V) is given by the righthand side of equa- 
tion ([h]). Here, I is viewed a parameter to be kept con- 
stant. A vector normal to E^ is then given by 

N A =d A *t(U,R,V). (23) 

If Ef is embedded in Minkowski space, then we find that 

N A N A = _4^(^-2)/(l-a) a (2-4a) c7 2a (~A £ ^ g M ^ 

(24) 



where we have used ([19]) to eliminate i? 2 , which means 
that N A N A is implicitly evaluated on E^. We note that 
N A is clearly null for U — 0, which corresponds to TVflE^, 
and spacelike for U > 0, which corresponds to the rest 
of E^. Hence, the components of the unit normal to E^, 
given by 



i A = N A l\N B Nu\ 1/ \ 



(25) 



diverge as U — » + . Since the 5D manifold is flat, the 
Riemann tensor on E^ is related to the extrinsic curvature 
4-tensor in the following manner: 



K a f3 = e A e B V B n A , 



R 



aft-yS 



(26) 
(27) 



where e A = dy A /dC, a are vectors spanning E^ and C" are 
the (arbitrary) coordinates used on E^. Because of the 
bad behaviour of n A at U = 0, we fully expect the extrin- 
sic curvature and the 4D Riemann tensor to be singular 



along TV n E^. Indeed, our extrinsic curvature formal- 
ism is not really defined at U = 0, so we must again 
be content with a limiting argument. We can calculate 
various curvature scalars composed from K a p (or equiv- 
alently R a f3-ys) and show that they diverge as U — » + . 
For example, consider 



K a ?K a0 = (V B n A )(V A n B ), 



(28) 



which can be established from the definition of K a $ and 



use of the induced metric h AB = h a pe A eg = g AB 
n A n B . For E^ 6 M5, we obtain via computer: 

(a 2 -2a + 2 ) a (2«-2)£(4a-2)/( Q -i) 



K afi K, 



a0 



u 2a 

as U -> 0+ (y A 6S/6 M 5 ).(29) 



Therefore, Ef has a curvature scalar singularity at all 
points with U = 0; i.e. all along the line TV n E^. This 
confirms the conclusion of the t — > + limiting argument 
used above. 

It is interesting to note that this curvature singularity 
is due to N A becoming null along TVfl E^, as opposed to 
any lack of "smoothness" of the surface at U = 0. To see 
this, we can consider the embedding of E^ in Euclidean 
space. Then, the norm of N A along TV fl E^ (i.e. [7 = 0) 
is 

N A N A = 2V 2 , [y A £ (TVn E £ ) e E 5 ]. (30) 

Unlike the Minkowski case, N A N A only vanishes at a 
single point on TV (~l E^. We also obtain via computer 



ae (0,i) 
a £ (|, 00) ' 



(31) 



where y A £ (TV n E^) £ E 5 . For a £ (|,oo), we get 
a curvature singularity for U = V = R = 0, which is 
the position of the caustic in Figures [l] - Interestingly 
enough, for a £ (0, |) we see no singularity at all, despite 
the fact that N A N A vanishes at U = V = R = 0. This 
is because the components of N A themselves vanish at 
U = V = R = 0, which means that the unit normal n A is 
finite and well defined (as may easily be verified). This is 
in agreement with Figures || and [|, which suggests that 
Ef is smooth at the origin. We have also investigated 
the behaviour of the full K a p 4-tensor evaluated along 
TV (~l E^ in E5 and confirmed that the components are in 
general well-behaved, except at V = for a £ (i, 00). 
This reinforces our conclusion that the only Euclidean 
curvature singularity onE^isat£/ = V r = i? = for 
a £ (4, 00). The chief difference between the Minkowski 
and Euclidean embeddings is that n A is infinite along 
TVnE^ in the former case, while in the latter scenario the 
unit normal is well defined all over E^. Therefore, we can 
attribute the Minkowski line-like curvature singularity on 
E^ to the divergence of the unit normal, and not to the 
lack of smoothness of the surface along TV D E^. 
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In conclusion, the geometric structure of the big bang 
in 5D depends on one's definition of a singularity. If a 
singular location on is taken to be a place beyond 
which timelike geodesies cannot be extended, then the 
big bang is a point-like event in 5D. If a singular location 
is taken to be a place where curvature scalars on di- 
verge, then the big bang manifests itself as a line in 5D. 
The distinction can be elucidated by asking whether an 
observer at r = at time t = to is ever in causal contact 
with points other than the 7 r caustic on TV PI . The 
trajectory of a light ray passing through r — at t = to 
is easily found from the null geodesic condition applied 
to the 4D metric (0) : 



-(*) = 



a ^- a /(l-a)]M-l)/o t _ £(a-l)/a] 

(a-1) 



(32) 



Substituting for r in (^) and taking the t — ► + limit, 
we find that light arriving at r = at any finite time to 
must have come from T = R = L = for a E (4, oo) and 
T = —L — > -co, R = for a G (0, i). In other words, 
the only point on the AfC\T,g line inside the cosmological 
horizon of r — is coincident with the position of the 
7 r caustic. This is true for all reasonable values of a 
and to. Because of isotropy, this must also be true for 
any finite value of r. Therefore, of all the points on 
where curvature scalars diverge, only the j r caustic is 
in causal contact with "normal" points on S^. In a real 
physical sense, this means that the geometric structure of 
the singularity that gives rise to the observable universe 
is point-like in 5D. Information from the rest of the line- 
like singularity on can never reach us within a finite 
amount of time. We feel that this is the best possible 
resolution to the issue of the geometric structure of the 
big bang in five dimensions. 



IV. THE ol — > AND a -> 1 LIMITS, AND THE 
a = \ PONCE DE LEON COSMOLOGY 

The coordinate transformation (g) between the origi- 
nal Ponce de Leon metric (|]) and the Minkowski metric 
( |l0| ) is undefined for a — 0, |, and 1. Of these three 
possibilities, only the a = | case leads to a non-singular 
5D metric tensor ([!]) in the x A coordinate system. We 
will first discuss the cases for which the Ponce de Leon 
cosmologies themselves are ill-defined and then turn to 
the highly-symmetric a = ^ case. 

In the limit a — > 0, equation (||) implies that the equa- 
tion of state of the induced matter is that of de Sitter 
space, namely p + p = 0. It is well known that the 
scale factor of the de Sitter universe grows exponentially 
in time, in contrast to the scale factor of the Ponce de 
Leon cosmologies which can only grow as fast as a power 
law. Therefore, it makes sense that de Sitter space cor- 
responds to the a — > limit of (|l|) since it is for this case 
that the scale factor oc t x l a grows faster than any (finite) 
power of t. On the other hand, the a — ► 1 case has an 



equation of state p + ip = 0, which represents matter 
with no gravitational density by equation (Q). The lim- 
iting value of the scale factor is oc t in this case, which 
is consistent with q(r) — » in equation (Q). This limit 
clearly corresponds to the empty Milne universe, which 
is known to have scale factor linear in time and matter 
with p + 3p = 0. It is interesting to note that q(r) = 
would correspond to the osculating paraboloid of be- 
ing cylindrical, which is precisely the situation interme- 
diate to the a = 3/2 case shown in Figure |l| and the 
a = 2/3 case shown in Figure |^. 

Now, as mentioned above, the a = ^ cosmology is the 
only case for which 5D metric (0) is well defined, but the 
transformation (|^) is not. This cosmology represents the 
case intermediate between universes with a 7 r caustic 
atT = i? = L = (Figures [I] - 0) and those with a 
7 r caustic at null infinity (Figures^ and |^). The line 
element by (0) is 



dS 2 =?dt 2 -t 4 £ 4 da 2 -t 2 d£ 2 



(33) 



By t — > it this gives a metric with signature ( h), 

so the 4D spacetime is of the type used in Euclidean 
quantum gravity [ pd| . A further change I — ► it gives 
back ((3^). The metric is also unchanged for t <-> I with 
one or the other of t — > it and I — > il. However, (|3^) with 
the last term taken with reversed sign is not a solution 
of the field equations Rab = 0, as may be verified by 
computer. It is therefore not like the so-called "two-time" 

5D metrics with signature (H h) which satisfy 

Rab = and have been applied to vacuum waves pQ] 
and cosmological inhomogeneities fiH . It can also be 
mentioned that the Ponce de Leon metric ([!]) does not 
satisfy Rab — if the last term is taken with reversed 
sign. Thus while (Q) and ( |33| ) satisfy the field equations 
with a space-like extra dimension, they do not with a 
time-like extra dimension, which raises the possibility of 
testing the signature of higher-dimensional theory. 

The model ( j33"l) also shares with ([!]) the property that 
the 5D Riemann-Chistoffel tensor is Rabcd = 0. This 
may be verified by computer. The scale factor in ( |3^ ) 
varies as t 2 , which is faster than standard FPW models. 
It contains 4D matter which by @ with r = it is a 
perfect fluid with density p — 3/2itt 2 and pressure p = 
— 1/ttt 2 . The equation of state is therfore p = — 2p/3. 
The gravitational density ([?]) is (p + 3p) = — p, so the 
model is inflationary. 



V. CONCLUSION 

iV-dimensional extensions of general relativity and par- 
ticle physics are now numerous and include Kaluza-Klein 
theory, induced-matter theory, superstrings, supergrav- 
ity, string theory and membrane theory. Among these 
possibilities, the Ponce de Leon models provide a con- 
crete realization of cosmology in 5D because they embed 
the FPW solutions of general relativity along with their 
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properties of matter. Mathematically, a smooth local em- 
bedding of 4D solutions with a finite energy-momentum 
tensor in a class of 5D solutions of the vacuum field equa- 
tions follows from Campbell's theorem. Physically, the 
4D solutions exist on hypersurfaces of the 5D manifold, 
and represent matter-filled spacetimes because the 4D 
space is curved, even though in the Ponce de Leon case 
the 5D space is flat. What we have done in the present 
account is to concentrate on the physical consequences of 
how curved 4D cosmologies are embedded in the flat 5D 
Ponce de Leon cosmologies. 

One of our main results is to use the algebraic trans- 
formation (|^) from the models in FRW coordinates ([!]) 
to Minkowski coordinates ( |l0| ) to obtain pictures of the 
big bang. Figures [l] - [H] elucidate the structure of sev- 
eral types of 4D universes by showing them as 2-surfaces 
embedded in a 3D space. The universes shown include 
standard models for late (matter-dominated) and early 
(radiation-dominated) epochs, as well as several infla- 
tionary cases. The congruence of comoving geodesies 
on any given E^ hypersurface is found to radiate from 
a caustic at t = Q. The caustic appears at different 5D 
positions for models with a £ (0, ^) and a € (^, oo), but 
it always falls on the null ray J\f. The curvature of the 
hypersurfaces diverges along the intersection of TV and 



S^, but observers on E^ are not in causal contact with 
any of the singular points in A/"nEf except for the caustic 
in the comoving congruence. This leads us to conclude 
that the big bang is point-like in 5D. The case (^) where 
our coordinate transformation is undefined has been in- 
vestigated numerically and is also inflationary. 

This and other aspects of our work would repay an 
in-depth study. The main restriction on the results pre- 
sented here is that they refer only to the Ponce de Leon 
solutions of the vacuum 5D field equations |2q| . There 
is an extensive literature on other solutions of the these 
field equations Q|. Some of these are not 5D flat like 
the Ponce de Leon cosmologies. However, we see no im- 
pediment to extending our approach using embeddings 
to curved 5D solutions and N(> 5)D solutions. 



Acknowledgments 

We thank A.D. Linde for discussions, NSERC for sup- 
port and CIPA for hospitality. We would also like to 
thank the anonymous referees, whose suggestions greatly 
improved the presentation of our results. 



S. W. Hawking and G. F. R. Ellis, The Large Scale 
Strucutre of Space-Time (Cambridge University Press, 
Cambridge, 1973). 

R. M. Wald, General Relativity (University of Chicago 
Press, Chicago, 1984). 

A. D. Linde, Inflation and Quantum Cosmology (Aca- 
demic Press, Boston, 1990). 

P. S. Wesson, Space-Time-Matter (World Scientific, Sin- 
gapore, 1999). 

A. D. Linde, Rep. Prog. Phys. 42, 389 (1979). 

R. Brout, Phys. Rev. Lett. 43, 417 (1979). 

A. H. Guth, Phys. Rev. D 23, 347 (1981). 

R. N. Henriksen, A. G. Emslie, and P. S. Wesson, Phys. 

Rev. D 27, 1219 (1983). 

W. B. Bonner, J. Math. Mech. 9, 439 (1960). 

P. S. Wesson, Astron. Astrophys. 151, 276 (1985). 

Euclidean Quantum Gravity, edited by G. W. Gibbons 

and S. W. Hawking (World Scientific, Singapore, 1993). 

A. Vilenkin, Phys. Lett. B 117, 25 (1982). 

F. Darabi, W. N. Sajko, and P. S. Wesson, Class. Quant. 

Grav. 17, 4357 (2000). 

J. E. Lidsey, C. Romero, R. Tavakol, and S. Rippl, Class. 
Quant. Grav. 14, 865 (1997). 

Modern Kaluza-Klem Theories, edited by T. Appelquist, 
A. Chodos, and P. G. O. Freund (Addison- Wesley, 
Menlo-Park, 1987). 

J. E. Campbell, A Course of Differential Geometry 
(Claredon Press, Oxford, 1926). 

S. Rippl, C. Romero, and R. Lavakol, Class. Quant. Grav. 
12, 2411 (1995). 

C. Romero, R. Tavakol, and R. Zalaletdinov, Gen. Rel. 
Grav. 28, 365 (1996). 



[19] E. Witten, Nucl. Phys. B 186, 412 (1981). 

[20] P. West, Introduction to Supersymmetry and Supergravity 

(World Scientific, Singapore, 1986). 
[21] M. B. Green, J. H. Schwarz, and E. Witten, Superstring 
Theory (Cambridge University Press, Cambridge, 1987). 
M. J. Du ff Int. J. Mod. P hys. A 11, 5623 (1996). 
D. Youm, [hep-th/0004144] (2000). 



[22 
[23] 

[24] R. Maartens Jiep-tIi/0004166 

[25] ' 
[26] 



(2000). 

A. Chamblin, fiep-th/001 1 L2g| (2000). 



H. Redmount, Mon. Not. 



D. Lynden-Bell, J. Katz, and 
R. astr. Soc. 239, 201 (1989). 
[27] W. Rindler, Phys. Lett. A 276, 52 (2000). 
[28] J. Ponce de Leon, Gen. Rel. Grav. 20, 539 (1988). 
[29] P. S. Wesson, Astrophys. J. 394, 19 (1992). 
[30] P. S. Wesson, Astrophys. J. 436, 547 (1994). 
[31] P. S. Wesson, Astrophys. J. 440, 1 (1995). 
[32] D. J. McManus, J. Math. Phys. 35, 4889 (1994). 
[33] A. A. Coley and D. J. McManus, J. Math. Phys. 36, 335 
(1995). 

[34] G. Abolghasem, A. A. Coley, and D. J. McManus, J. 

Math. Phys. 37, 361 (1996). 
[35] K. Lake, P. J. Musgrave, and D. Pollney, GR Tensor 

(Department of Physics, Queen's University, Kingston, 

Canada, 1995). 

[36] H. Liu and B. Mashhoon, Ann. Phys. (Leipzig) 4, 565 
(1995). 

[37] P. S. Wesson, H. Liu, and S. S. Seahra, Astron. Astro- 
phys. 358, 425 (2000). 

[38] M. Lipshutz, Theory and Problems of Differential Geom- 
etry, Schaum's Outline Series (McGraw-Hill, New York, 
1969). 

[39] A. D. Linde, Scientific American 271, 48 (1994). 



[40] A. Billyard and P. S. Wesson, Gen. Rel. Grav. 28, 129 
(1996). 

[41] A. Billyard and P. S. Wesson, Phys. Rev. D 53, 731 



